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ABSTRACT 

We apply the Union2 compilation of 557 supernova la data, the baryon acoustic oscillation 
measurement of distance, the cosmic microwave background radiation data from the seven 
year Wilkinson Microwave Anisotropy Probe, and the Hubble parameter data to study the ge- 
ometry of the Universe and the property of dark energy by using models and parametrizations 
with different high redshift behaviours of w(z). We find that ACDM model is consistent with 
current data, that the Dvali-Gabadadze-Porrati model is excluded by the data at more than 3a 
level, that the Universe is almost flat, and that the current data is unable to distinguish mod- 
els with different behaviours of w(z) at high redshift. We also add the growth factor data to 
constrain the growth index of Dvali-Gabadadze-Porrati model and find that it is more than la 
away from its theoretical value. 
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1 INTRODUCTION 

Even since the discovery of the accelerated expansion of the Uni- 
verse in 1998 dRiess et al.ll 19981 : iPerlmutter et alj|l999h , many ef- 
forts have been made to confirm and understand this phenomenon 
of acceleration. For the explanation of the acceleration, there are 
three different approaches. The first method introduces a new ex- 
otic form of matter with negative pressure, dubbed as dark energy 
to drive the Universe to accelerate. The cosmological constant is 
the simplest candidate of dark energy which is also consistent with 
observations, but at odds with quantum field theory. The second 
method modifies the theory of gravity known as general relativ- 
ity at the cosmol ogical scale, such as the Dvali-G abadadze-Porrati 
(DGP) model dDvali. Gabadadze & Porratill2000h . The third ap- 
proach takes the view that the Universe is inhomogeneous. In this 
paper, we focus on dark energy and DGP models. 

The only effect of dark energy we know is through gravita- 
tional interaction; this makes it difficult to understand the phys- 
ical nature of dark energy. In particular, the question whether 
dark energy is the cosmological constant remains unanswered. 
Recently, it was claimed that the flat ACDM model is inconsis- 
tent with observations at more t han l a level (Huang et al 120091 : 



ing behaviour of dark energy. However, no evidence for dark en- 
ergy dynamics was found in other s t udies (Lampeitl et al. 2009 
Serraetalj|2009l : iGong et alj[2010al : iGong. Wang & Cail |20 10b 



2009; ICai, Su & Zhand l20ld) . 



Shafieloo, Sahni & Starobinsky 

Shafie loo. Sahni & Starobinskvl ([2009) suggested that the cosmic 
accele ration is slowing down from z ~ 0.3. In iHuang et al.l 
(2009), it was c laimed that dark energy suddenly emerged at red- 
shift z ~ 0.3. ICai. Su & Zhand J2010h found possible oscillat- 
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Pan et al. 2010). The diffe rence betwee n the conclusion s drawn 

in lShafieloo. Sahni & Starobinskvl fe009l) and lGong etai] ( l2010ij) 
lies in th e Baryon Acoustic Oscillation (BAO) data used in their 
analysis. Shafieloo, Sahni & Starobinsky ( 2009) employs the ra- 
tio Dy (0.35) / Dy (0.2) of the effe ctive distance D v {z) at two 
redshifts, while Gong et al. (2010a) applies the BAO A parame- 
ter given by lEisenstein et ajj J2005I) . The tension between BAO 
measureme nt and higher r edshif t type la supernova (SN la) was 
no ticed in |P ercival et al. d2007h . and the tension was lessened 
in IPercival et alj d20ld) due to revised error analysis, different 
methodology adopted and more data. 

It was argued that the sy stematics in different data sets heav- 
ily affected the fitting results {H icken et al. 2009; Sollerma n et al.l 
l2009l:lGong. Wang & Caill2010bl : lKessler et al.ll2010h . The Consti- 
tution compilation found that the scatter at high redshift is higher 
for SALT and SALT 2 fitters, and SALT 2 poorly fits the nearby 
U-band light curves faickenetal.ll2009h . However, it was found 
that SALT2 performs better than SALT and MLCS2k2 judged 
by the scatter around the b est-fi tting luminosity distan ce relation- 
ship in lConlev et alj d2008l) and [Amanulla h'etal] d201dh - Because 
MLCS2k2 training is mainly based on the observation of nearby 
SN la, and because the observations made in the observer-frame 
U-band are contaminated with a high level of uncertainty due to 
atmospheric variations, so MLCS2k2 is less accu rate at predicting 
the rest-frame U-band usi ng high redshift SN la lAmanullah et alj 
l20ld : iKessler et al.ll2010u . The Union2 data applies the SALT2 
light curve fitter because it addresses the problem by including 
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high redshift data where the rest-frame U -band is observed at red- 
der wavelengths dAmanuilah et alj|20ld) . In this paper, we com- 
bine the Union2 sampl e of 557 SN la data with systematic errors 
dAmanuilah et al.ll201oh . the BAO distance ratios r s (z d )/Dv(z) 
between the comoving sound horizon at the baryon drag epoch 
r s (zd ) and the effective distance Di/(z) at z — 0.2 and z = 
0.35 dPercival et ai]|20ld). the radial BAO measurements at z 



0.24 and z = 0.43 dGaztanaga. Miguel & Sanchez! l2009ah . the 



seven-year Wilkins on Mi crowave Anisotropy Probe (WMAP7) 
data l Komatsu et al. 201 lh. and the Hubble param eter H(z) data 



dGaztanaga. Cabre & Huil2009bl ; IStern et al]201Cl) to probe the ge- 
ometry of the Universe and the nature of dark energy with different 
models. 

The paper is organis ed as follows. I n sec tion 2, we 
present the SN la data dAmanuilah et al J l20ld). the BAO 
data dGaztanaga. Miguel & Sanchez! 2009al; [percival et al.l 



the WMAP7 data ( Komatsu et al. 



the 



data (Gaztanaga, Cabre & Huil l2009bl: ISternetalJ l2010h. and 


the growth factor data (Viel, Haehnelt & Springel 2004 




2006; 


McDonald et ajj 20051; iTegmark el al] 120061: Ross et al 




2007; 


Angela et al. 2008; Guzzo et al] 2008: Blake et al. 20ld). 


and all 



the formulae related to these data. In section 3, we list all the 
models and the fitting results, and conclusions are given in section 
4. 



2 OBSERVATIONAL DATA 

The Union2 SN la data consist of the low-z SN la data observed 
at the F.L. Whipple o bservatory of the Harvard-Smithsonian cen- 
tre for astrophysics dHicken et al J 120091) . the intermediate-z data 
observed during the first seas on of the Sloan Di gital Sky Survey 
(SDSS)-II supernova survey dKessler et al. 201 0J), and the high- 
z data from the Union compilation (Kowalski et al. 2008). The 
Union2 SN la data used the SALT2 light-curve fitter because it 
performs better than both SALT and MLCS2k2 when judged by 
the s catter around the bes t-fitting luminosity distance relation- 
ship dAmanuilah et al.ll2010h . To use the 557 Union2 SN la data 
dAmanuilah et alfcOlOl) . ~ 



we minimize 



= ^ M Zi ) ~ IAobs{Zij\C s n{Zi,Zj)[v{Zi) - )J-obs{zj 



(1) 



where the extinction-corrected distance modulus /i(z) = 
51og 10 [dL(z)/Mpc] + 25, C sn (zi,Zj) is the covariant ma- 
trix which includes th e systematical errors for the SN la data 
dAmanuilah et al.ll2010h ; the covariant matrix is available onlinfj 
The luminosity distance cIl(z) is 



d L (z) 



l + z 



■■ Sk 



(2) 



the dimensionless Hubble parameter E(z) — H(z)/Ho; and 
Sk{x) is defined as x, sin(a;) or sinh(:r) for k = 0, +1, or -1, re- 
spectively. Due to the arbitrary normalization of the luminosity dis- 
tance, the nuisance parameter h = _ffo/(100 km s _1 Mpc _1 ) in the 
SN la data is not the observed Hubble constant. So we marginalize 
the nuisance parameter h with a flat prior, after the marginalization, 
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we get dGong. Wu & Wandl2008l) 

XL(P) = o<iC~J(zi, Zj)oy 




(3) 



where a* = fJ, a bs(zi) — 25 — 51og 10 [I/odz,(z;)]- 

In addition to the Union2 SN la data, we use the BAO distance 
measurement from the oscillations in the distribution of galaxies. 
The BAO is due to the sound waves in the plasma of the early Uni- 
verse and the wavelength of the BAO is related to the comoving 
sound horizon at the baryon drag epoch. The distance at the red- 
shift z — 0.2 was measured in the clustering of the combined 2dF 
Galaxy Redshift Survey (2dFGRS) and SDSS main galaxy sam- 
ples, and the distance at the redshift z = 0.35 was measured in 
the clustering of the SDSS luminous re d galaxies. From the BAO 
observation of the galaxy power spectra. fpercival et al.l (2010) mea- 
sured the distance ratio, 



d 2 = 



r s {zd) 
D v (z) 



(4) 



at two redshifts z = 0.2 and z = 0.35 to be dg I = 0.1905 ± 
0.0061, and d%% = 0.1097 ± 0.0036 (hereafter Bao2). Here the 
effective distance is 



D v (z) = 



d 2 L (z 



(l + z) 2 H(z) 



1/3 



(5) 



Zd is the drag redshift defined in lEisenstein & Hul dl998h , the co- 
moving sound horizon is 



r s (z) 



c s (x)dx 



(6) 



where the sound speed c s (z) = 
R b = 3Q. b h 2 /{4 x 2.469 x 10~ 
calculate 



l/y/3[l + R b /(1 + z)], and 
To use this BAO data, we 



1 ^CzUdud^Adj 



XBao2(P,V b h ,h) = 



(7) 



where di — (d z= o.2, ^2=0.35), Ac^ — di — d° bs and the covari- 
ance matrix CBao(di, di) for the two paramete rs do. 2 and do. 35 is 
taken from equation (5) in lPercival et al.l d2010f) . Besides the model 
parameters p, we need to add two more nuisance parameters Q,bh 2 
and Q,mh 2 when we use the BAO data. 

From the measurement of the radial (line-of-sight) BAO scale 
in the galaxy power spectra, the cosmological parameters may be 
determined from the measured values of 



AzBao{z) 



H(z)r e (z d ) 



(8) 



at two redshifts z = 0.24 and z = 0.43, which are 
Az B ao{z = 0.24) = 0.0407 ± 0.0011 and Az Bao (z = 
0.43) = 0.0442 ± 0.0015 (her eafter Baoz), respectively 
(Gaztanaga, Miquel & Sanchez 2009a). Therefore, we add x 2 with 

XBaozKPAlbH ,tl) - { {hOOTl ) ( o\ 

. f Az Bao (0.43)-0.0442 \ 2 ( ' 

+ V 0.0015 J ■ 

When we add this BAO data to the fitting, we also need to add 
the nuisance parameters f2j,/i 2 and Q. m h 2 . In iGong. Wang & Cail 
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d2010bh . it was found that the Baoz data is consistent with the Bao2 
data, and the constraints on the model parameters get improved 
with the addition of the Baoz data. 

Because both the SN la and the BAO data measure the distance 
up to redshit z < 2, we need to add distance data at higher redshift 
z > 10 to better constrain the property of dark energy, so we use 
the WMAP7 data. When the full WMAP7 data are applied, we need 
to add some more parameters which depend on inflationary models, 
and this will limit our ability to constrain dark energy models. So 
we only use the WMAP7 measurements of the derived quantities 
such as the shift parameter R(z*) and the acoustic scale Ia(z*) at 
the decoupling redshift, and the decoupling redshift z* . Then we 
add the following term to "x 2 , 



2 

XCMB 



{p,Q b h 2 ,h) = ^ ^-CBt^Ai)^ 



(10) 



where the three parameters Xi = [R(z*), Ia(z*), z*], Axi = 
Xi — x° and the covariance matrix CcMBjxj, Xj) for the three 
parameters is taken from Table 10 in Komatsu et alj d201 lh . The 
shift parameter R is expressed as 




R{z*) 



The acoustic scale I a is 



ttcIl(z*) 



(1 + z*)r a {z*) 



dz 

eH) 



302.09 ± 0.76, 



1.725±0.018.(11) 



(12) 



an d z* is the decoupling r edshift with the parametrization defined 
in iHu & Sugivamal dl996l) . We also need to add the parameters 
Qbh 2 and Q m h 2 to the parameter space when we employ the 
WMAP7 data. 

The SN la, BAO and WMAP7 data measured the distance 
which depends on the double integration of the equation of state 
parameter w(z), the process of double integration smoothes out 
the variation of equation of state parameter w(z) of dark energy. 
To alleviate the double integration, we also apply the measure- 
ments of the Hubble parameter H(z) which depends onQ.DE di- 
rectly and detects the variation of w(z) better than the distance 
scales. Furthermore, the add ition of the H(z) d ata can better con- 
strain w(z) at high redshift dGong et al]|2010j) . In this paper, we 
use the H(z) data at 11 different redshift s obtained f r om t he dif- 
ferential ages of red-envelope galaxies in lStern et alj yOlO), and 
three more Hubble parameter data H(z = 0.24) = 76.69 ± 2.32, 
H(z = 0.34) = 83.8 ± 2.96 and H(z = 0.43) = 86.45 ± 3.27, 
determined by iGaztanaga. Cabre & Huil J 2009b). So we add these 
H(z) data to % 2 , 



(13) 



where o^i is the lcr uncertainty in the H(z) data. Basically, The 
model parameters p are determined by minimizing 



2 _ 2 , 2 

X Xsn ' XBaol 



2 . 2 

XBaoz + XCMB 



2 

XH- 



(14) 



In addition, we add the prior h = 0.742±0.036 dRiess et al. 2009). 

In order to distinguish the modified gravity such as DGP 
model from dark energy models, we approximate the growth factor 
with/(z) = Q 1 m + { 1 -A/7)Q k taong. Ishak & Wand2009h . then 
we use the growth factor data f(z) obtained from the measurement 
of t he redshift distortion to constrain t he growth index 7 of the mod- 
els dViel. Haehnelt & Springelll2004l2006l : lMcDonald et"ai]|2005l ; 



i Tegmark elakl 120061; IRoss et al.1 120071; lAngela etafl 120081; iGond 
l2008l : lGuzzo et al.ll2008l ; lBlake et al.l20ld : lDossett et al.|2010h . So 

we calculate 



X/(P>7) 



15 

E 

i = l 



[/(%)-n^(^)-( 7 -4/7)n fc ( Zi )] 2 



(15) 



7i 



The likelihood for the parameters p in the model and the 
nuisance parameters is computed using the Monte Carlo Markov 
Chain (MCMC) method. The MCMC method randomly chooses 
values for the above parameters p, evaluates \ 2 and determines 
whether to accept or reject the set of parameters p using the 
Metropolis-Hastings algorithm. The set of parameters that are ac- 
cepted to the chain forms a new starting point for the next process, 
and the process is repeated for a sufficient number of steps until 
the required convergence is reached. Our M CMC code is based on 
the publicly available pa ckage COSMOMC dlxwis & Bridlell2002l : 
iGong, Wu & Wang||2008l) . 

After fitting the observational data to different 
dark energy models, we apply the Om diagnostic 

to detect the de- 
For 



dSahni. Shafieloo & Starobinskvl l2008h 
viation from the ACDM mode l. 
dSahni. Shafieloo & Starobinskvll2008h . 

Om(z)= E ^- X 



a flat universe 



(16) 



(l + 2) 3 -l' 

For the ACDM model, Om(z) = Q m is a constant which is in- 
dependent of the value of f2 m . Because of this property, Om diag- 
nostic is less sensitive to observational errors than the equation of 
state parameter w(z) does. On the other hand, the bigger the value 
of Orn(z), the bigger the value of w(z). 



3 COSMOLOGICAL MODELS 
3.1 qi — <?2 parametrization 

To understand the current accelerating expansion, we parametrize 
the de celeration paramete r q(z) with a simple two-parameter func- 
tion dGong & Wang 2003), 



, , 1 qiZ + q? 

g(2) =2 + (TT7F 



(17) 



In this parametrization, we have only two parameters p = (qi , q^). 
The parameter q2 = q(z = 0) — 1/2, and q(z) = 1/2 at high red- 
shift which represents the matter dominated epoch. In principle, 
this parametrization does not involve Q. m and Qk, but the comov- 
ing distance depends on the geometry of the Universe through the 
function Sk, for simplicity, we consider the flat case fifc = for 
this model. Although the flat assumption of Qk = may induce 
large error in the estimation of cosm ological parameters due to 
the de generacy among O m , D.k and w dClarkson. Cortes & Bassettl 
2007), but for this model, the only effect of Qk is through Sk, and 
Sk(x) w x when Qk is small, so the impact of the flat assumption 
is small. The dimensionless Hubble parameter is 



E(z) 



exp[/ z [l + gH]dln(l + u)| 



(1 + z) 3/2 exp 



+ 



27T+7T 2 " 



(18) 



When z > 1, E 2 (z) w (1 + z) 3 exp(gi + q 2 ), so we can think 
171+92 = lnf2 m . To account for the radiation-dominated Uni- 
verse, we take the above E(z) as the approximation for the matter 
and dark energy only, so we use the following Hubble parameter to 
approximate the whole expansion history of the Universe, 
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E\z) 



Q r {l + zf + (1 + zfexp 



qi + 



qiz 



<12 



(1 + -) 5 



(19) 



where the current rad iation component Q r — 4.1736 x 10~ 5 /i~ 2 
(Kom atsu et alJ[2oTTh . Fitting the model to the combined SN la, 
Bao2, Baoz, WMAP7 and H(z) data, we get the marginalized la 
constraints, qi = 0.07 ± 0.11 and q 2 = -1.43 ± 0.09 with x 2 = 
542.6. In terms of qo = q(z — 0), we find that go < —0.62 
at 3cr confidence level, so the evidence of current acceleration is 
very strong. Furthermore, we find that the 3a constraint on Q m is 



0.257+ 



The contour plot for Q m and qo is shown in 



Fig.CJd). 

At a low redshift, the radiation is negligible, so Om(z) in this 
model is 

Qm(z) Jl + ^^ + (^-^)/(l^) 2 ]-l (20) 
w (1 + zf - 1 

By using the best-fitting values of gi and q^, we reconstruct Om(z) 
and the results are plotted in Fig.[2jd). 



3.2 Piecewise parametrization of q(z) 

To further study the evolution of the deceleration parameter q(z), 
we use the more model-independent piecewise parametrizations. 
We group the data into four bins so that the number of SN la in 
each bin times the width of each bin is around 30, i.e., NAz ~ 30 
and N = 4. The four bins are zi — 0.1, 22 = 0.4, 23 = 0.7, 
24 = 1.8 and 25 extends beyond 1089. For the redshift in the range 
Zi-i ^ 2 < Zi, the deceleration parameter q(z) is a constant qi, 
q(z) — qi. Take 20 = 0, then for 2;_i < 2 < Zi, we get 



E(z) = (! + , 



(21) 



In this model, we have four parameters p = (qi, qz, qs, qi) - In 
general, for the piecewise parametrisation, the parameters such as 
qi in different bi ns are correlated and their errors depend upon each 
other. We follow lHuterer & Cooravl {2005) to transform the covari- 
ance matrix of g^s to decorrelate the error estimate. Explicitly, the 
transformation is 



Qi = yVijgj, 



(22) 



3.3 ACDM model 

For the cosmological constant, the equation of state parameter w = 
p/p = —1, and the energy density pa is a constant. For a curved 
ACDM model, the curvature term Qk 7^ 0, Friedmann equation is 

E 2 (z) 



fi fc (l + zf + Q m (l + zf + 0.(1 + zf + «A. (23) 



At low redshits, the contribution of the radiation term is negligible. 
We have two parameters p = (Q m , Qk) and one nuisance param- 
eter h in this model. By fitting the ACDM model to the combined 
SN la, Bao2, Baoz, WMAP7 and H(z) data, we get the marginal- 
ized ltr constraints, Q m = 0.272±° and n ft = 0-002 ± 0-004 
with \ 2 = 541.2. The contours of Q m and Qk are plotted in 
Fig. EJa). By fitting the model to observational data combined 
with the growth factor data, the marginalized la constraints are, 
Q m = 0.272±°;of , Q h = 0.002 ± 0.004 and 7 = 0.56± ;£| with 
X 2 = 546.3. 



3.4 DGP model 

In this model, gravity appears four di mensional at short distances 
while modified at large distances jDvali. Gabadadze & PorratH 
2000). The model is motivated by brane cosmology in which our 
Universe is a three-brane embedded in a five dimensional space- 
time. The Friedmann equation is modified as 



E 2 (z) = Q k {l + zf + p d + ^fi m (l + 2) a + fi 



212 



(24) 



where Q<j = (1 — Q m — f2fc)/2-\/l — If we take the point of 
view that Friedmann equation is not modified and the extra term 
in equation d24t is dark energy, then the equivalent dark energy 
equation of state parameter 111(2) for the DGP model is 



w{z) 



041 + zf + 2Q d [y / fl m {l + zf + Q 2 d + fl d 



2[On(i + zf + + n dy /Q m {i + zf + n 2 d ] 



(25) 



when 2 > 1, w(z) 1/2 and w(z = 0) = -(1 - Q*,)/(l + 

$7 m — $7fc). Since 0^ is very small, w(z) > — 1 for the DGP model. 

By fitting the DGP model to the combined SN la, Bao2, 
Baoz, WMAP7 and H(z) data, the marginalized la constraints are 
n m = 0.288l^; lf and Q k = 0.019 ± 0.005 with x 2 = 561.6. 
By fitting the DGP model to the combined SN la, Bao2, Baoz, 
WMAP7, H(z), and f(z) data, we get the marginalized la estima- 



tions Q„ 



0.290 



+0.014 



n fe = 0.019±0.005, and 7 



withx = 567.5. 



46+ 12 



where the transformation matrix T = V T r~ 1 ^ 2 V, the orthogonal 
matrix V diagonalizes the covariance matrix C of qi and T is the 
diagonalized matrix of C. For a given i, Ty can be thought of as 
weights for each gj in the transformation from qi to Qi. We are free 
to rescale each Qi without changing the diagonality of the correla- 
tion matrix, so we then multiply both sides of the equation above 
by an amount such that the sum of the weights . T%j is equal to 
1. This allows for easy interpretation of the weights as a kind of 
discretized window function. Now the transformation matrix ele- 
ment is Tij I Tit and the covariance matrix of the uncorrelated 
parameters is not the identity matrix. The i-th diagonal matrix el- 
ement becomes Tif)~ 2 . In other words, the error of the un- 
correlated parameters Qi is <7, = 1/ ^\ Tij . Fitting the model to 
the combined SN la, Bao2, Baoz, WMAP7 and H(z) data, we get 
the constraints on the uncorrelated parameters Qi and the result is 
shown in Fig. [3] 



3.5 CPL parametrization 

For the Chevallier-Polar s ki-Linder (C PL) parametrization 
dChevaliier & Polarskill200ll : lLindedl2003l) . the equation of state 
parameter is 



w(z) — Wo + ■ 



(26) 



so 10(2 = 0) — wo and w(z) ~ wo + w a when 2 3> 1. The 
corresponding dimensionless dark energy density is 

n DE {z) = 0,(1 + z ) 3 ( 1 +»0-H"a) e l-3™ a s/{l+*)] ! (27) 

where fl x = 1 — Q m — £l r — Qk- In this model, we have four 
model parameters p = (fl m , Qk, wo, w a ). Fitting the model 
to the combined SN la, Bao2, Baoz, WMAP7 and H(z) data, we 
get the marginalized la constraints, Q. 

(-0.26 
-0.06 



0.008±° ;°°s, wo 



1.16l^, and w a 



- 0-265^Q Qog, r^fc 
0.691?;^ with x 2 
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540.5. The contours of wo and w a are plotted in Fig.[TJa), and the 
contours of Q m and Ofe are plotted in Fig.[4jb). 
For the flat CPL model, Om(z) becomes 



Om(z) 



ft m (i + z) 3 + M^)-i 



(28) 



(l + z) 3 -l 

where Qde(z) is defined in equation {27} with Qk = 0. By fit- 
ting the combined data to the flat CPL model, we get the marginal- 
ized la constraints, Q, m = 0.267±%%{ 9 , w = -lmto] 7 , and 
w a = 0.07to'.88 with x 2 = 541.1. Using this result, we recon- 
struct Om(z) with equation d28t and the result is shown in Fig. 
Ha). 



3.6 JBP parametrization 

For the Jassal-Bagla-Padm a nabha n 
jjassal. Bagla & Padmanabhanl |2005|) . 
parameter is 



j(z) = w + 



WqZ 

W+7) 



(JBP) parametrization 
the equation of state 



(29) 



so w(z — 0) = wo and w(z) ~ wo when z 3> 1. In this model, 
the parameter wo determines the property of the equation of state 
parameter w(z) at both low and high redshifts. The corresponding 
dimensionless dark energy density is then 



3(l+u>o) /2(l+z) 2 ] 



(30) 



where O x = 1 — Q. m — Q r — fij.. In this model, we also have 
four parameters p = (f2 m , flk, wo, w a )- Fitting the model to 
the combined SN la, Bao2, Baoz, WMAP7 and H(z) data, we 
get the marginalized la constraints, Q m = 0.263l?,'n 2 , fijb = 

H0. 32 



0.004 ± 0.006, wo = -1.21] 



and w a = 1-29 



o.oi 

+ 1.35 



with 



X 2 = 540.6. The contours of Q, m and flfc are plotted in Fig.[4]c), 
and the contours of wo and w a are plotted in Fig.[]2b). 
For the flat JBP model, Om(z) becomes 



Om(z) 



n m (l + z) 3 + Q DE (z) 
(1 + z) 3 - 1 



(31) 



where Qde(z) is defined in equation J30t with Q, k ~ 0. By fit- 
ting the combined data to the flat JBP model, we get the marginal- 
ized la constraints, Q. m = Q.265±%° {1, w = — 1.08±g;?|, and 
w a — 0.32^i'72 with \ 2 = 541.0. Using this result, we recon- 
struct Om(z) with equation ( Bit and the result is shown in Fig. 

Hb). 
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Figure 1. The marginalized la and 2a contour plots of wo (f2 m ) and w a 
(go) for the CPL (a), JBP (b), Wetterich (c) and q(z) (d) parametrisations. 
The dashed line in the upper left panel denotes the condition Wo + w a = 0. 



are plotted in Fig.[TJc), and the contours of Q. m and fife are plotted 
in Fig.Hd). 

For the flat Wetterich model, Om(z) becomes 

Om { z)=^ 1+ f: + ^~\ (34) 



(1-M) 3 -1 

where Q.de (z) is defined in equation i33\ with Q k = 0. By fitting 
the combined data to the flat Wetterich model, we get the marginal- 
ized la constraints, Q. m = 0.266± o ; o i 5 , w o = -1.05l o ;?g, and 
w a = 0.14 ± 0.1 with x 2 = 541.1. Using this result, we recon- 
struct Om(z) with equation l !34t and the result is shown in Fig. 

He). 



3.8 Piecewise parametrization otw(z) 

Finally, we consider a more model-independent parametrization of 
w(z), the piecewise parametrization of w(z). In this parametriza- 
tion, the equation of state parameter is a constant, w(z) = Wi for 
the redshift in the range Zi-i < z < Zi, For convenience, we 
choose 20 = 0. We also assume that w(z > 1.8) = — 1. For a flat 
Universe, if Zi-i ^ z < 



3.7 Wetterich parametrization 

Now we consider the parametrization proposed by IWettericr] 
(t2004h . 



w(z) 



Wo 



(32) 



[l + w a ln(l + 2)] 2 ' 

For this model, w(z = 0) = wo and w(z) ~ when z > 1, so 
the behaviour of w(z) at high redshift is limited. The dark energy 
density is 

V DE =(l-O m -Q k - Or)(l + z )8+3«*/P-H..In(l+.)] < (33) 

In this model, the model parameters are p = (fi m , Q k , Wo, w a )- 
Fitting the model to the combined SN la, Bao2, Baoz, WMAP7 
and H(z) data, we get the marginalized la constraints, Q m = 
0.264 ± 0.013, Q k = 0.009+ 014 



0.32t°;« with x 2 



o.oos^ ««o = -1.17± :£. and 
540.4. The contours of too and w a 



n D E{z) = (i-n m )(i+z) 



Hd+z 



\3(u>i- 



H) 



(35) 



Again, the four parame ters Wi are correlated and we follow 
Huter er & Cooravf d2005l) to transform these parameters to decor- 
related parameters Wi. By fitting the model to the combined SN 
la, Bao2, Baoz, WMAP7 and H(z) data, we get the error estima- 
tions of Wi and the results are shown in Fig. [5] 



4 CONCLUSIONS 

We summarize all the results in Table 1 and some results are shown 
in Figs. lllSl By parametrizing the deceleration parameter q(z), we 
find very strong evidence for the current acceleration. For the piece- 
wise parametrization of q(z), we find that q(z) < in the redshift 
range ^ z < 0.6, and q(z) > in the redshift range z > 0.8 as 
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Table 1. The marginalized lcr errors for parameters constrained by the observational data 





X 2 /DOF 




n k 


WO (<?l) 


Wa (<?2> 


AIC 


BIC 


ACDM 


541.2/576 


n 979 +0.013 
-0.011 


0.002 ± 0.004 






545.2 


553.9 


DGP 


561.6/576 




0.019 ± 0.005 






565.6 


574.3 


CPL 


540.5/574 


265+ 019 


0.008,°;^ 


-1 16+ - 26 


0.69±?; 2 « 


548.5 


565.9 


JBP 


540.6/574 


0.263^0^ 


0.004 ± 0.006 


i 91 +0.32 
x - z -0.18 


i.29ij : |§ 


548.6 


566.0 


Wetterich 


540.4/574 


0.264 ±0.013 


009+ - 014 
u.uu»_ 005 


i 17 +0.09 
11 ' -0.23 


u, ° -0.16 


548.4 


565.8 


Qi — 92 model 


542.6/576 






0.07 ±0.11 


-1.43 ±0.09 


546.6 


555.3 


flat CPL 


541.1/575 


267+ ' 019 
u - ZD '-o.oi 




-i-oslo 1 ;! 7 


u - u '-0.88 


547.1 


560.2 


flat JBP 


541.0/575 


265+ ' 019 




-l-OBtss 


0.32+^ 


547.0 


560.1 


flat Wetterich 


541.1/575 


266+ 01 
u.^oo_ 015 






0.14 ±0.1 


547.1 


560.2 




Figure 3. The lcr and 2cr errors of the four uncorrelated Qi. The red dashed 
line is reconstructed with the best-fitting ACDM model. 



Figure 5. The lcr and 2cr estimates of the four uncorrelated parameters Wi . 
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shown in Fig.[3] So the Universe experiences accelerated expansion 
up to the redshift z ~ 0.6 and decelerated expansion at large red- 
shift z > 0.8. For the CPL, JBP and Wetterich models, we see from 
Fig.Q]that ACDM model is consistent with them, and this is further 
confirmed by the Om diagnostic as shown in Fig. [2] The piecewise 
parametrization of w(z) also confirms that ACDM model is con- 
sistent with current observations as shown in Fig. [5] The CPL, JBP 
and Wetterich models differ in the behaviour of w(z) at high red- 
shift; from Table 1 we see that all of them fit the observational data 
well, and w(z) < as seen in Fig.QJa). So the current data is still 
unable to distinguish models with different behaviours of w(z) at a 
high redshift. The number of parameters for ACDM and DGP mod- 
els are the same, the difference between the best-fitting value of \ 2 
is A% 2 = 20.4, so DGP model is excluded by the current data at 
more than 3<r level. The observational constraint on the growth in- 
dex 7 is 7 = 0.46lo.J|g for tne DGP model which 

is more than 

la away from the theoretical value 11/16, and the growth index 
of ACDM model is 7 — 0.56^q'qq which is consistent with the 
theoretical value 6/11. Our results also show that the Universe is 
almost flat. By using the prior —5 ^ log |flfc| ^ 0, it was found 
that -0.9 x 10 -2 sC tt k sC 0.01 at 9 9% confidence level with 
the B ayesian model averaging method ( Vardan van. Trotta & SiUd 
1201 lh . In order to compare different models with different num- 
ber of parameters, w e usually apply Akaike Information Criterion 
(AIC) i lAkaikd[l974T) . In terms of AIC, ACDM model is slightly 
preferred by the current data. Furthermore, to account for the ef- 
fects of the number of data points and the number of parameters, 
Bayesian Information Criterion (BIC) dSchwad[l974h is used for 
model comparison. In terms of BIC, the ACDM model is again pre- 
ferred by the current data. In addition to the approximate methods 
like AIC or BIC for model comparison, the Baye sian mode l com- 
parison provides a better tool for model selection dTrottj2008l) . 

Our results are based on the stan dar d \ 2 method whic h has 
some shortcomings dMarch et alj201 lh . so lMarch et alj 1201 ll) pre- 
sented the Bayesian hierarchical method to constrain the cosmo- 
logical parameters and argued that the new method gives tighter 
constraint and outperforms the standard \ 2 method. 
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